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Abstract 


\ 

/ 

The  technique  of  pole  placement  was  applied  to  a  time  varying 
periodic  system.  When  the  control  matrix  was  of  full  rank,  unstable 
real  roots  were  predictably  moved  to  stable  locations.  Specified 
damping  was  successfully  applied  to  oscillatory  roots.  In  the  case 
of  the  control  matrix  having  less  than  full  rank  a  least  squares  gain 
was  applied.  This  minimized  the  error  between  all  desired  pole  loca¬ 
tions  and  controlled  pole  locations  and  as  a  consequence,  prediction 
of  final  pole  location  was  not  possible.  An  exact  method  for  pre¬ 
dictably  placing  poles  as  desired  was  restricted  by  a  controllability 

condition.  i  '4. 

/ 

The  method  was  applied  to  the  restricted  four-body  problem  and 

/ 

control  of  a  satellite  in  a  periodic  orbit  about  the  Earth-Moon 
Lagrangian  point,  L^,  A  stabilized  orbit  was  obtained  for  a  model 
which  allowed  instantaneous  control  of  both  position  and  velocity. 
Stabilization  was  not  successful  in  the  case  of  control  of  velocity 
only.  The  least  squares  solution  effected  several  modes  simultaneously. 
The  exact  solution  for  gain  was  unobtainable  as  the  controllability 


condition  was  not  met. 


I.  INTRODUCTION 


Background 

The  phenomenon  of  gravitational  interactions  between  satellites  and 
celestial  bodies  has  long  been  a  point  of  interest.  In  the  past  several 
decades  the  understanding  of  this  phenomenon  has  become  a  requirement  for 
satellite  control.  In  nature,  interactions  between  multiple  large  body 
systems  can  generate  points  of  collected  small  objects.  An  example  is 
the  Trojan  objects  in  the  Sun-Jupiter  system.  The  phenomenon  is  the 
physical  realization  of  the  periodic  n-body  problem  in  dynamics. 

The  Trojan  objects  were  predicted  by  the  solutions  of  the  restrict¬ 
ed  three-body  problem  as  put  forth  by  Lagrange  (1772).  In  this  problem, 
three  bodies  exhibit  coplanar  motion.  One  has  mass  so  small  its  motion 
does  not  effect  the  motion  of  the  remaining  two;  this  motion  being  cir¬ 
cular  orbits  rotating  at  a  constant  rate  about  the  barycenter.  Lagrange 
found  five  equilibrium  solutions  to  the  equations  of  motion.  These  were 
points  where  gravitational  forces  on  the  small  third  body,  as  exerted 
by  the  other  two,  would  balance  the  "centrifugal"  force,  thus  leaving  it 
fixed  with  respect  to  a  reference  frame  rotating  with  the  two  large 
bodies.  Figure  I  presents  the  rotating  frame  for  the  Earth-Moon  system. 

A  body  put  at  one  of  the  Libration  points  with  zero  relative  velocity 
would  stay  there  until  perturbed.  Gylden  called  the  Lagrange  solutions 
"centers  of  Libration"  since  in  this  rotating  frame  the  third  body  may 
also  describe  small  time-periodic  orbits 


about  the  fixed  points.  The  equilibrium  positions  have  been  designated 
Lj  thru  Lg  and  are  shown  in  figure  2  for  the  Earth-Moon  system.  The 
rotating  frame  has  its  origin  at  the  Earth-Moon  center  of  mass.  The 
points  L1  to  L3  are  col  linear  with  the  Earth-Moon  line  of  sight  while 
and  Lg  each  make  an  equilateral  triangle  with  the  Earth  and  Moon 
forming  two  vertices.  Motion  near  the  col  linear  points  has  been  shown 
to  be  unstable  through  infinitesimal  analysis  (Ref  5:420-428).  Unstable 
motion  infers  that  a  body  at  these  points  would  depart  if  perturbed. 

The  triangular  points  show  infinitesimally  stable  motion  in  their  neigh¬ 
borhood  from  a  similar  analysis,  meaning  the  bodies  at  and  Lg  would 
remain  in  that  vicinity  even  when  perturbed.  Discovery  of  the  Trojan 
objects  which  appear  to  librate  near  the  Sun-Jupiter  Lagrange  points 
gives  impetus  to  this  result. 

In  the  literature,  distinction  is  made  between  stability  at  the 
libratior.  point  itself  and  stability  of  the  periodic  solutions  about 
the  point,  if  they  exist  (Ref  12:215-217).  The  fundamental  questions 
at  hand  about  motion  near  the  Lagrange  points  concern  the  existence 
of  periodic  orbits  and  their  stability.  For  the  restricted  three-body 
problem,  these  questions  have  been  answered  (Ref  5,8,12,14).  Un¬ 
fortunately,  in  implementing  an  actual  satellite  the  use  of  the 
restricted  three-body  problem  solutions  in  a  control  scheme  does  not 
suffice.  Addition  of  a  fourth  body  alone  without  effects  of  eccen¬ 
tricity  or  orbital  inclination  render  the  triangular  points  unstable. 
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Figure  2  Lagrange  Points  in  Earth-Moon  System 


The  unstable  orbit  is  made  more  unstable  in  this  restricted  four- 
body  problem. 

Militarily,  the  stability  and  control  of  periodic  systems, 
specifically  the  libration  point  orbits  has  moved  from  the  realm  of 
academic  interest  to  a  requirement  for  the  continued  secure  use  of 
space  for  surveillance,  early  warning,  command,  communication,  and 
control.  Current  satellites  for  those  purposes  are  located  in  low 
and  near  Earth  orbits.  The  Soviet  Union  has  repeatedly  demonstrated 
its  ability  to  intercept  and  destroy  low  Earth  orbit  satellites. 

Recent  revelations  have  indicated  the  Soviets  now  possess  the  capa¬ 
bility  to  attack  even  satellites  in  geosynchronous  orbit  using  an  ICBM 
as  a  direct  ascent  attack  vehicle  (Ref  11).  The  move  to  the  libration 
points  orbits  would  have  not  only  the  advantage  of  increasing  intercept 
time-of-fl ight  from  hours  to  days  but  has  the  advantage  of  increasing 
the  area  of  coverage  of  each  satellite.  Synchronous  satellites  in 
orbit  about  the  Earth's  equator  are  at  an  altitude  such  that  the 
North  Pole  is  over  the  horizon  and  hidden  from  view.  Thus,  this  is  a 
blind  spot  in  a  surveillance  net  of  geosynchronous  satellites  requir¬ 
ing  many  low  orbit  satellites  with  high  orbital  inclinations  for 
continuous  coverage  instead.  However,  the  libration  point  orbits  are 
in  the  Earth-Moon  orbital  plane  and  hence  at  any  time  one  satellite 
can  have  a  direct  line  of  sight  to  the  Pole.  Also,  with  satellites 
located  at  the  libration  points,  global  communications  would  be  possibl 
through  one  or  more  relays  to  these  satellites.  Thus,  we  have  a  vested 
interest  to  understand  and  solve  the  problem  of  control  of  periodic 


systems. 

Problem  and  Scope 

The  thrust  of  this  effort  is  the  control  of  unstable  modes  of  a 
periodic  system.  A  controllability  statement  on  oscillatory  time 
varying  modes  and  unstable  real  modes  will  be  attempted.  By  this  we 
want  to  state  under  what  conditions  can  we  "stabilize"  the  system  and 
arrange  system  poles  predictably.  To  do  this,  an  attempt  will  be  made 
to  adapt  pole  placement  techniques  for  linear  constant-coefficient 
systems  to  linearized  periodic  systems.  The  use  of  singular  value 
decomposition  along  with  techniques  of  generalized  inverses  will  be 
explored  to  find  the  appropriate  system  gain  required  to  stabilize  the 
system.  The  restricted  four-body  problem  will  be  used  as  a  practical 
application  for  the  techniques  developed. 

Previous  work  has  centered  on  efficiently  controlling  a  satellite 
affected  by  the  unstable  real  mode  of  the  l_3  periodic  orbit  as  developed 
by  Wiesel  (Ref  16).  The  length  of  time  control  has  been  effective  at 
reasonable  cost  has  been  increased  for  the  restricted  four-body  problem 
but  still  is  not  long  enough  for  practical  application.  Smith  (Ref  10) 
applied  linear,  constant  gain,  feedback  control  to  a  satellite  near  the 
orbit  calculated  by  Wiesel.  The  control  law  he  developed  was  inadequate 
based  upon  the  failure  of  various  gains  and  feedback  quantities  used  to 
return  the  satellite  to  the  reference  orbit  economically.  Shelton's 
(Ref  9)  use  of  modal  control  allowed  satellite  control  for  up  to  two  years 


II.  FLOQUET  THEORY 

The  Hamiltonian  for  the  four-body  problem  devised  in  reference  9 
results  in  equations  of  motion  which  are  explicit  functions  of  time. 

The  variational  equations  of  this  nonautonomous  system  can  be  written 
in  linearized  state-variable  form 

6x  =  [A(t)]6x  (1) 

The  elements  of  the  nxn  matrix  A(t)  are  continuous  periodic  functions 
of  time  with  period  T  such  that  the  periodicity  is  expressed  as 

C A ( t ) ]  =  [A( t+T ) ]  (2) 

As  a  consequence  of  this  periodicity,  we  will  be  able  to  predict 
motion  for  all  time  once  the  solution  for  system  motion  is  found  over 
one  period.  To  do  this  requires  using  Floquet's  Theorem.  This  theorem 
concerns  the  representation  of  a  fundamental  system  of  solutions  of  (1) 
as  the  product  of  a  periodic  matrix  (with  period  T)  and  a  solution 
matrix  for  a  system  with  constant  coefficients.  The  fundamental  solu¬ 
tion  matrix  has  been  given  the  symbol  4>(t)  and  can  be  expressed  as 

t[R] 


O(t)]  =  [Q(t)]e 


(3) 


Here,  Q(t)  is  periodic  where 


CQ(t)]  =  [Q( t+T ) ]  (4) 

and  the  columns  of  i<(t)  are  related  to  the  columns  of  Q(t)  by  e*^;  R 
a  constant  matrix.  The  fundamental  matrix  satisfies 

[+(t)3  =  [A(t)]  [*(t)]  (5) 


If  *l> ( t+T )  satisfies 


[[’(t+T)]  =  [A( t+T ) 3  [<t>(t+T)]  ,  (6) 

then  by  (2),  it  also  satisfies 

[♦(t+T)]  =  [A ( t ) ]  [ 4> ( t+T ) ]  (7) 

and  is  therefore  itself  a  fundamental  matrix.  Since  both  4>(t)  and 
•i>( t+T )  are  solution  matrices  containing  n-linearly  independent  solu¬ 
tions  to  (1),  the  columns  of  the  matrices  must  be  related  such  that 


[♦(t+T)]  =  [*(t)]  [C] 


(8) 


where  the  matrix  C  is  called  the  monodromy  matrix.  For  t=0,  equation 
(8)  can  be  rearranged  to 


O-1(0)]  O(T)]  =  [C]  (9) 

which  will  be  used  later.  The  monodromy  matrix  is  not  unique  but 
depends  upon  the  particular  fundamental  matrix  used  to  satisfy  (5).  A 
matrix  <i>(t)  which  satisfies  (5)  but  also  is  the  identity  matrix  at 

H 

t^O  is  called  the  principal  fundamental  matrix.  The  principal  funda¬ 
mental  matrix  gives 


0(T)p]  =  [C]  (10) 

using  (9)  and  from  the  fact  that  the  inverse  of  an  identity  matrix  is 
itself. 

We  now  evaluate  equation  (3)  at  t=0  and  take  the  inverse  of  the 
resulting  matrix  equation.  Combining  the  outcome  of  this  with  equa¬ 
tion  (3),  evaluated  now  at  t=T,  we  get 

[+■'(0)]  O(T)]  =  [Q_1(0)]  [Q(T)]eT[R:l  (11) 

Since  Q(t)  satisfies  (4)  the  product  [Q"*(0)]  [Q(T)]  reduces  to  the 


identity  matrix  and  by  invoking  (9),  equation  (11)  is  reduced  to 


L2 


The  eigenvalues  of  the  C  matrix  are  called  the  characteristic 
multipliers  of  A(t)  and  are  symbolized  by  A..  (i=l,2, . . .  ,n) .  The 
characteristic  multipliers  are  related  to  the  R  matrix  eigenvalues,  p. 


Pi  =  f  n  (A.)  (14 

The  p.  are  alternately  called  the  characteristic  or  Poincare  exponents 
The  Poincare  exponents  indicate  the  system  stability  with  positive 
real  (or  real  part  for  complex)  eigenvalues  being  unstable.  Since  any 
monodromy  matrix  for  a  particular  problem  is  similar  to  all  others 
(Ref  5:266)  and  hence  has  the  same  eigenvalues,  the  characteristic 
multipliers  as  well  as  the  characteristic  exponents  can  be  put  in 
Jordan  canonical  form.  When  all  the  eigenvalues  of  a  matrix  are 
unique,  the  Jordan  form  is 


In  the  case  where  some  of  the  eigenvalues  are  complex  conjugate  pairs, 
a  modified  block  diagonal  form  of  the  Jordan  matrix  is  used.  This 
allows  the  creation  of  an  eigenvector  corresponding  to  the  real  part 
of  the  oscillatory  mode  and  an  eigenvector  for  the  imaginary  part, 
replacing  the  two  complex  eigenvectors  of  the  conjugate  roots.  The 
block  diagonal  [J]  matrix  is  of  the  form 


[J] 


- 

°1  "  “i  •  •  •  0 

U>J  •  0 

o  o  x3  .  . 


0  0  .  0  An 


(16) 


where  a  is  the  real  part  and  u>  the  imaginary  part  of  the  complex  root. 

There  are  both  real  and  complex  Poincare  exponents  for  the  period¬ 
ic  orbit  about  1^  (see  Table  I)  so  (16)  depicts  the  Jordan  form  of  R. 
Using  this  in  (3)  and  (5)  we  get 


[<3(t)]et[J]  +  [Q(t)J  =  [A(t)J  [Q(t)]et[J]  (17) 
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TABLE  I 


Normalized  Poincare  Exponents  of  Shelton's  Orbit 


Poincare 


Exponent 


Planar  Mode 

0.0 

+ 

0.077738i 

0.0 

- 

0.077738i 

Planar  Mode 

0.17893 

+ 

O.Oi 

-0.17893 

+ 

O.Oi 

Out  of  Plane 

0.0 

+ 

0.084134i 

0.0 

_ 

0.084134i 

The  e  L  J  terms  are  removed  leaving 

COU)]  +  [Q(t)]  [J]  =  [A(t)]  [Q(t)] 
which  can  be  solved  for  [J] 

[J]  =  [Q“l(i-)]{tA(t)]  [Q(t)]  -  [Q(t)  ]} 
Now,  defining  A{t)  as  Q(t)  partitioned  into  columns 

[  A(t)]  =  [(Ut)  , 


and  putting  this  into  (19)  gives 

[0]  -  [  Al(t)]  {[A(t)3  «t)]  -  [Wt)]J  (21) 

The  matrices  a  and  a” ^  are  periodic  with  period,  T.  Remembering 
the  linearization  of  the  variational  equations  we  had 


ox(t)  =  x(t)  -  Xq( t ) 


3x(t) 

ax(tQ) 


x0(t) 


(  tg  )  J 


(22) 


stating  the  variations  were  the  difference  between  the  actual  and  a 
reference  solution.  The  reference  solution  used  here  being  the  Weisel 
periodic  orbit.  Defining 


[<fr(t,tQ)]  = 


3x(t) 

Dx(tQ) 


XqU) 


we  get  (22)  into  the  form 


(23) 


«x(t)  =  C4>(t,t0)]6x(t0)  (24) 

where  ,<» ( t , t0 )  is  the  state  transition  matrix  having  the  properties 
[+(t,t0)]  =  [♦(t,t1)]  [♦(t1,t0)] 

(25) 

[■t'(t0,t0)]  =  [I]. 
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Since  i(t)  are  solutions  to  (1)  as  is  6x(t),  (24)  becomes 


<f4 


[v(t)]  =  [>t(t,t0)]  [4*(t0)l  ,  (26) 


which,  for  ty=0  and  t=T  we  get 


[+(T)3  =  [4>(T,0)]  [$(0)]  (27) 


solving  for  ^ ( T ,0 )  and  combining  with  equation  (9)  we  can  show 


[C]  =  |>(T,0)] 


(28) 


An  arbitrary  solution  to  (1),  which  is  a  column  from  (3),  can  be 
written  as 


6x.(t)  =  Q^tJePi1  (29) 

Equation  (24)  at  tQ=0  and  (29)  equate  such  that 


6x(t)  =  Q.^^i4  =  [♦(t,0)]6xi(0).  (30) 


Now,  with  t-0  in  (29) 


Ax^t)  =  Qi(t^pit  =  [♦(t,0)]Q.(0)  (31) 
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and  again  using  (4) 


P,T  PiT 

■’ii-d)  =  Qi(T)e  1  =Qi(0)e  =  ♦(T,0)Qi(0)  (32) 

at  the  end  of  one  period. 

Equation  (32)  is  reduced  to  an  eigenvalue  eigenvector  problem 
using  (28) 

{[C]  -  e"]t  [I]]  Qi(0)  =  0  (33) 

,,.t 

where  the  terms  e  are  the  eigenvalues  and  Q^(0)  the  eigenvector.  Tc 
get  a  complete  solution  the  eigenvector  must  be  integrated  over  an 
entire  period.  This  is  done  by  getting  Aon  the  left-hand  side  of  (21) 
We  can  then  solve  (21)  numerically  since  it  is  really  an  ordinary  dif¬ 
ferential  equation  in  l{ t)  where  these  are  just  the  partitioned  Q(*’ 
and  its  initial  conditions  are  the  eigenvectors  of  C.  The  values  of 
A-*(t)  are  obtained  directly  from  A(t)  and  A(t)  is  composed  of  n-lin- 
early  independent  vectors,  is  of  full  rank  and  is  invertible.  Again, 
since  the  A  and  /"*  matrices  are  periodic,  once  the  solution  is 
obtained  over  one  period  it  is  known  for  all  time. 


III.  POLE  PLACEMENT 


Pole  placement  in  linear  constant  coefficient  systems  allows  the 

designer  to  specify  system  time  response  through  the  relocation  of 

specific  system  modes  or  poles.  Alteration  of  the  time  response  is 

limited  by  system  zeroes.  This  is  accomplished  by  selective  application 

of  gain  in  the  control  law.  The  technique  leaves  non-specified  modes 

unchanged.  The  technique  can  be  applied  to  any  controllable  mode,  hence 

the  requirement  for  a  control labil ity  statement.  Also,  the  system  can 

be  stabilized  only  if  uncontrollable  modes  are  stable.  Stabilization 

is  accomplished  by  moving  poles  from  the  right-half  of  the  complex 

plane  (unstable)  to  the  left-half  plane.  Damping  may  be  added  to  a 

complex  pair  of  roots  (again  moving  the  roots  to  the  left)  but  this 

must  be  accomplished  by  moving  pairs  of  complex  conjugates  (Ref  3:362). 

* 

For  Hamiltonian  systems,  the  Poincare  exponents  are  always 
positive/negative  pairs.  This  tells  us  that  if  we  have  a  real  root  in 
the  left-half  plane  its  conjugate  is  in  the  right-half  plane  and  the 
system  is  unstable.  Complex  roots  also  appear  as  conjugates.  The 
problem  at  hand  has  two  oscillatory  modes  and  a  real  mode.  The  real 
mode  lies  in  the  orbital  plane  and  has  one  unstable  root.  The  oscil¬ 
latory  modes,  one  mode  in  and  one  mode  out  of  the  plane,  are  purely 
imaginary  (marginally  stable)  in  the  restricted  four-body  problem  but 
become  unstable  in  more  accurate  models. 


The  technique  of  pole  placement  will  be  used  to  move  the  unstable 
real  root  to  a  specific  point  in  the  left-half  plane  along  with  adding 
damping  to  the  oscillatory  modes. 

Implementation  of  pole  placement  for  control  requires  augmenting 
the  state  equation  such  that 


6x(t)  =  [A(t)]Sx(t)  +  [B(t)]u (t) 


Here,  the  matrix  B  applies  control  to  the  desired  states  and  u  is  the 
control  vector.  The  system  requires  transformation  to  a  new,  uncoupled, 
set  of  coordinates  to  insure  any  control  applied  has  only  the  desired 
effect,  specifically,  leaving  the  stable  real  pole  unaltered.  This  new 
set  of  modal  variables  is  generated  by 


n  =  [  A_1(t)]6x(t). 


The  state  vector  is  transformed  by  solving  (35)  for  Sx(t)  and  taking 
its  time  dei ivative,  yielding 


6x(t)  =  [ A(t)]n(t)  +  [A(t)]n(t) 


letting  (34)  equal  (36)  and  again  using  the  solution  of  (35)  for  6x(t) 
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we  have 


ff» 


C  A(t)].|(t)  +  [  A(t)]ri(t)  =  [A ( t ) ]  [  A(t)]n(t)  +  [B(t)]u(t)  (37) 


which  is  readily  solved  for  7,(t) 


'i 


(t) 


t /_1  (t)  J  {[A(t)]  [  A( t ) ]  -  [ *A(t)]]  n(t)  +  [  a”  1  ( t ) 3  [B(t)]u(t) 

(38) 


Remembering  equation  (18)  and  equation  (20) 


I'*  t)3  ^  [A(t)]  [  a(  t )  ]  -  [  A(  t )  ]  [J]  (39) 


and  by  defining 


5(0  =  CG(t)3  C a( t ) ]  n(t) 


(40) 


we  derive 


~'(t)  =  [J]  7,(t)  +  [  A_1(t)]  [B(t)]  [  A(t)]n(t)  (41) 

Upon  rearranging,  this  becomes 

n(t)  =  {[J]  +  [  A~ 1  ( t ) ]  [ B ( t ) ]  [G( t ) ]  [  A(t)]j  7,(t)  (42) 
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(43) 


M  *  [J]  +  [  4-1(t)]  [B (t)]  [G(t)3  ( A ( t ) 3 

we  write 

7,(t)  =  [K]n(t)  (44) 

The  crux  of  the  problem  is  we  have  picked  K,  the  location  of  system 
poles  and  must  now  solve  for  the  gain  matrix,  G,  that  gives  us  this 
desired  response.  This  is  done  by  first  obtaining 

[B]  [G]  =  [A(t)]{[K]  -  [J]}  [  A_1(t)]  (45) 

then  solving  for  the  gain. 

The  matrices  K,  J,  and  A  are  all  known.  The  difficulty  arises  in 
finding  the  solution  to  (45)  when  B  is  singular  or  is  a  nonsquare  matrix. 
The  first  problem  we  will  address  is  the  one  where  we  have  a  square  B 
matrix  of  full  rank.  Since  this  matrix  is  of  full  rank  it  is  non¬ 
singular  and  invertible.  The  solution  for  the  gain  matrix  can  be  found 
directly  from 

[G(t)]  =  [B_1(t)]  [A(t)]  [[K]  -  [J])  [  A- 1  ( t ) ]  (46) 

The  system  in  this  case  is  always  controllable  since  the  product  [B] 

[G]  in  equation  (45)  has  no  zero  row  (Ref  1:444). 


Two  other  methods  for  finding  the  gain  matrix  for  nonsquare  B 
matrices  were  used.  These  vary  in  the  way  a  solution  to  equation  (45) 
is  found.  An  optimum,  but  not  exact,  solution  is  obtained  using  the 
generalized  inverse  of  B  in  equation  (46)  for  the  standard  inverse.  In 
this  problem,  any  of  the  methods  of  obtaining  the  generalized  inverse 
in  the  next  section  may  be  used.  An  exact  solution  can  be  obtained 
using  the  method  of  singular  value  decomposition  and  its  resulting 
factorization  of  B  into  square,  invertible  matrices  and  a  matrix  of 
singular  values. 

The  optimum  solution  is  limited  in  that  it  results  in  the  construc¬ 
tion  of  a  least  squares  "best"  gain  matrix.  This  gain  does  not  predict¬ 
ably  move  a  specific  system  pole  to  the  desired  location.  It  moves  all 
system  polos  in  such  a  way  as  to  minimize  the  error  between  the  desired 
and  actual  pole  location.  The  result  is  that  the  system  may  or  may  not 
be  stabi 1 ized. 

Tiie  exact  solution  is  highly  desireable  but  the  application  of 
this  technique  is  limited  since  a  controllability  condition  must  be  met 
(see  Section  V).  As  a  result,  the  exact  solution  method  may  not  control 
the  system. 
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IV.  GENERALIZED  INVERSE 

Moore-Penrose  Pseudoinverse 

A  method  for  finding  the  generalized  inverse  of  an  mxn  matrix  B 
was  discovered  by  the  mathematicians  Moore  and  Penrose  earlier  this 
century.  Their  simple  formula  found  the  optimum  solution  of  the  matrix 
equation 


[A]x  =  b  (47) 

for  the  case  when  A  was  not  invertible  by  standard  means,  specifically 
when  A  was  not  square.  In  other  words  this  method  allows  us  to  "in¬ 
vert"  a  matrix  that  had  rank,  r,  equal  to  its  smallest  dimension, 
either  m  or  n. 

The  Moore-Penrose  formula  for  the  generalized  inverse  is 

B+  =  (BTB)  V  (48) 

T  T  _1 

where  B  is  the  standard  transpose  of  B  and  (B  B)  is  the  inverse  of 

the  nxn  square  matrix  of  full  rank  found  by  performing  the  multipli¬ 
cation  indicated.  When  B  itself  is  square,  B+  =  B*. 

Different  formulations  were  developed  for  B  having  a  greater 
number  of  rows,  m,  than  columns,  n,  (m>n)  and  for  the  case  of  m<n. 
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The  set  up  of  the  problem  in  this  paper  results  in  m  always  being 
greater  than  n,  so  only  this  case  will  be  discussed.  References  6  and 
13  detail  the  procedures  for  the  generalized  inverse  when  m<n. 

A  factorization  of  B  which  allows  us  to  create  a  generalized  in¬ 
verse  for  the  case  when  the  rank  of  B  is  less  than  the  smaller  of  m  or 
n  has  been  developed.  This  variation  of  the  Moore-Penrose  generalized 
is  termed  the  pseudoinverse  by  Strang  (Ref  13)  and  is  used  for  examples 
where  B^B  is  singular  and  therefore  not  invertible. 

B+  =  U,T(U'U,T)  Vl)  1L,T  ,  (49) 

and  a  proof  is  found  in  Strang  (Ref  13). 

The  terms  L'  and  U‘  are  from  the  factorization  of  B 

B  =  L'U* ,  (50) 

and  the  superscripts  T  and  -1  denote  the  standard  matrix  transpose  and 
inverse,  respectively.  The  matrices  B,  L',  U*  are  sequentially,  mxn, 
mxr,  rxn  with  all  three  having  rank,  r.  Matrix  U  is  formed  first  by 
Gauss  reduction  of  B  to  upper  triangular  form.  When  the  reduction  is 
done,  U  has  rank  r  with  the  bottom  m-r  rows  being  zero.  We  form  U*  by 


simply  eliminating  the  zero  rows  of  U.  Gauss  elimination  is  accomplished 
by  performing  elementary  row  operations  on  B.  This  is  actually  done  by 
premultiplying  elementary  matrices  onto  B  (Ref  13:20-21).  The  matrix 
L  is  formed  by  multiplying  the  inverses  of  the  elementary  matrices  in 
the  reverse  order  that  the  row  operations  were  performed.  We  reduce  L 
to  L'  by  eliminating  the  last  m-r  columns  of  L.  This  is  possible  since 
only  the  last  m-r  rows  of  U  would  multiply  these  and  they  are  zero. 
Singular  Value  Decomposition 

The  requirement  to  invert  non-square  real  matrices  has  lead  to  the 
development  of  several  techniques  now  in  wide  use.  One  of  the  lesser 
used,  though,  is  singular  value  decomposition.  An  important  point  of 
this  technique  is  that  it  allows  us  to  invert  a  rectangular  matrix  of 
less  than  full  rank.  A  matrix  such  as  this  would  occur  if  we  decided 
to  not  apply  control  to  one  mode  leaving  a  row  of  zeroes  in  B.  Like¬ 
wise,  as  with  the  pseudoinverse,  methods  have  been  developed  for  sys¬ 
tems  in  which  B  has  a  greater  number  of  rows  than  columns  and  also  for 
the  case  where  mn.  Again,  this  paper  always  has  m  greater  than  n  so 
the  reader  is  referred  to  reference  7  when  dealing  with  the  case  of  m<n. 

In  general,  any  mxn  matrix  B,  with  m>n,  can  be  factored  in  the 

form 

[B]  Mo]  [RET]  =  [LE]ID1[RET] 


(51) 


and  the  generalized  inverse  is  defined  as 


It  is  constructed  using  the  following  details.  The  matrix  LE  is  an 
mxm  orthogonal  matrix  and  its  columns  are  the  eigenvectors  of  BB^. 
Similarly,  the  columns  of  RE  are  the  eigenvectors  of  B^B  and  RE  is  nxn. 
Both  RE  and  LE  are  full  rank.  The  term  o  is  composed  of  two  submatrices, 
the  diagonal  D  matrix  containing  the  m  singular  values  and  the  remaining 
submatrix  is  mx(n-m)  with  zero  elements.  If  the  rank  of  the  matrix  is  r 
then  there  are  m-r  singular  values  that  are  zero.  The  singular  values 
are  the  positive  square  roots  of  the  eigenvalues  of  BBT,  that  is, 

4.  =  +[ci(BBT)]1/2  ,  i  =  1,2, ...,m  (53) 

therefore, 

D-diag  (Cj,  C2»*..tp»  0,...0)  (54) 

In  constructing  the  generalized  inverse,  B+,  the  only  step  of 
equation  (52)  needing  clarification  is  D+.  Since  D  is  diagonal  its 
inverse  is  simply 


D+  =  diag  (1/Cj.  l/c2*  —  .l/^r)- 


The  generalized  inverse  satisfies  the  following  four  conditions 


with  B 


(B+B)  =  B4B 


(BB+)  BB+ 


BB  B  =  B 


B+BB+  =  B+ 


Specific  properties  of  the  generalized  inverse  are  found  in  references 
6,  7,  and  13.  The  use  of  the  singular  value  decomposition  towards  this 
problem  along  with  numerical  examples  are  given  later. 
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.  CONTROLLABILITY  CONDITION 


In  the  section  on  pole  placement,  controllability  was  discussed 


briefly.  We  are  interested  in  solving  (45) 


[B]  [G]  =  U(t)]  j[K]  -  [J]j  [A_1(t)] 


for  the  gain  matrix  which  will  allow  us  to  predictably  relocate  system 
poles.  The  controllability  statement  is  the  restriction  on  the  extent 
to  which  we  can  predict  pole  location  and  place  poles  as  desired. 

When  B  has  full  rank  the  solution  to  (45)  is  straightforward  and 
G  is  found  at  any  time.  The  controllability  is  that  the  product  [B]  [G] 
have  no  zero  rows. 

For  the  case  where  B  is  nonsquare  with  m>n,  the  method  of  singular 
value  decomposition  can  lead  to  an  exact  expression  for  the  gain  provid¬ 
ed  a  consistency  check  (controllability  condition)  is  met.  In  the 
previous  chapter,  B  was  factored  in  (51)  to 


[B]  =  [LE]  [RE1]. 

0 

This  result  is  substituted  into  equation  (45)  giving 


[LE]  [RET]  [G(t»  =  [A(t)]([K]  -  [J]j  [  A_1(t)3 


Since  the  matrix  LE  is  of  full  rank  it  is  invertible.  It  is  also  orth¬ 


onormal  and  its  inverse  is  equal  to  its  transpose.  Premultiplying  both 
sides  of  (57)  by  LE  transpose  results  with 


D 

0 


[RET]  [G ( t ) ]  =  [LET]  [  a{ t ) ]  f[K]  -  [J]]  [  A- 1  ( t ) ]  =  C F ( t ) ] 


(58) 


The  matrices  F  and  G  can  be  written  in  vector  form 


[F(t)]  =  [F(t)1#  ?(t)2,...,  f(t)m]  (59) 


and 


[G(t)]  =  [g(t)p  g(t)2 . 9(t)J  (60) 


Equating  the  columns  of  (58) 


D 

0 


[Ret]  g(t)x 


7(t)i ,  i=l,m 


(61) 


We  can  solve  for  the  columns  of  the  gain  matrix  from 


m-r  , 

g  (t)  =  l  -1  (LE..f(t).)RE.  (62) 

1  j  =  l  r'j  3  J 


which  is  from  Noble  (Ref  6:338).  This  solution  is  found  provided 


This  i 


LC-f.  =  0  for  j=m-r, . . . ,m 
j  i 

and  i=l,m. 


s  also  the  condition  for  when  the  system  is  controllable. 


RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 


The  overall  scheme  for  implementing  the  controller  was  identical 
for  the  full  rank  B  matrix,  least  squares  method  for  nonsquare  B  matrix, 
and  exact  singular  valued  decomposition  of  B  techniques  in  determining 
the  gain  matrix.  Once  the  gain  matrix  was  obtained  by  one  of  these 
methods,  equations  (40)  and  (35)  were  used  to  generate  equation(34), 
which  was  then  solved  numerically.  The  Poincare  exponents  were  then 
obtained  and  system  stability  and  pole  location  checked. 

The  first  set  of  cases  examined  in  the  stabilization  of  a  satellite 
in  periodic  orbit  about  L3  was  for  B  having  full  rank.  Here,  B  was  a 
six-by-six  identity  matrix  allowing  control  to  be  applied  to  all  state 
variables  equally.  The  gain  matrix  as  a  function  of  time  was  calculat¬ 
ed  from  (46)  for  a  selection  of  pole  locations.  The  values  of  desired 
pole  location  as  well  as  the  resultant  Poincare  exponents  are  tabulated 
in  Table  II.  Values  for  desired  pole  location  were  picked  as  represen¬ 
tative  of  what  could  typically  be  chosen.  In  the  three  cases  done,  pole 
location  was  exactly  as  desired.  As  predicted,  with  all  modes  control¬ 
lable,  placement  of  poles  was  straightforward. 

This  result  is  important  for  a  periodic  system  where  all  the  states 
are  controllable.  The  technique  results  in  stable  control  in  a  pre¬ 
dictable  manner.  In  the  physical  terms  for  the  restricted  four-body 


problem,  however,  this  case  cannot  be  implemented  as  it  requires  instant 
aneous  control  over  satellite  position  (teleportation)  and  we  can  only 
control  momentum. 

TABLE  II 

Normalized  Poincare  Exponents,  Desired  Location 
and  Resultant  Location  for  Square,  Full  Rank  B  Matrix 


Case  Desired  Location  Resultant  Location 


Move  unstable 

0.0 

+  0.077738i 

0.0 

+  0.077738i 

real  pole 

0.0 

-  0.077738i 

0.0 

-  0.077738i 

-0.1 

+  O.Oi 

-0.1 

+  O.Oi 

-0.17893 

+  O.Oi 

-0.17893 

+  O.Oi 

0.0 

+  0.084134i 

0.0 

+  0.084134i 

0.0 

-  0.084134i 

0.0 

-  0. 084134 i 

Move  unstable 

-0.1 

+  0.077738i 

-0.1 

+  0.077738i 

real  pole. 

-0.1 

-  0.077738i 

-0.1 

-  0.077738i 

damp  first 

-0.1 

+  O.Oi 

-0.1 

+  O.Oi 

oscillatory 

-0.17893 

+  O.Oi 

-0.17893 

+  O.Oi 

mode 

0.0 

+  0.084134i 

0.0 

+  0.084134i 

0.0 

-  0. 084134 i 

0.0 

-  0. 084134 i 

Case 

Move  unstable 
real  pole, 
damp  both 
oscillatory 
modes 


TABLE  II  (Cont'd.) 


Desired  Location 
-0.1  +  0.077738i 

-0.1  -  0.077738i 

-0.1  +  O.Oi 

-0.17893  +  O.Oi 

-0.1  +  0  -  084134 i 


Resultant  Location 
-0.1  +  0 . 0777381 

-0.1  -  0.077738i 

-0.1  +  O.Oi 

-0.17893  +  O.Oi 

-0.1  +  0.084134i 


-0.1 


-  0.084134i 


-0.1 


-  0.084134i 


Use  of  this  result  directly  in  (45)  for  8"*  gave  the  time  varying 

gain  matrix  which,  when  implemented  into  the  system,  minimized  the 

error  between  desired  and  resultant  pole  location.  The  case  considered 

was  the  repositioning  of  the  unstable  real  root  to  the  same  location  as 

* 

was  used  for  the  full  rank  B  system.  The  normalized  Poincare  exponents 
desired  and  those  resulting  from  this  method  are  listed  in  Table  III. 


TABLE  III 

Normalized  Poincare  Exponents,  Desired  Location  and 
Resultant  Location  for  Nonsquare  B  Matrix  Using 
Generalized  Inverse  Method  of  Gain  Calculation 

Desired  Location  Resultant  Location 


0.0 

+ 

0.077738i 

-0.0845 

+ 

0. 1204i 

0.0 

- 

0.077738i 

-0.0845 

- 

0. 1204 i 

-0.1 

+ 

O.Oi 

0.0492 

+ 

0. 1979i 

-0.17893 

+ 

O.Oi 

0.0492 

- 

0.1979i 

0.0 

+ 

0. 0841 34 i 

0.0 

+ 

0.084134i 

0.0 


0.084 1341 


0.0 


0.084134i 


The  desired  relocation  was  not  obtained.  Both  the  real  mode  and  first 


oscillatory  mode  were  assigned  new  values  such  that  the  system  is  still 
unstable.  The  form  of  the  gain  matrix  and  coupling  between  planar  modes 
resulted  in  unpredictable  resultant  pole  location.  The  second  oscilla¬ 
tory  mode  was  unchanged  since  the  gain  matrix  had  no  corresponding  ele¬ 
ments  to  effect  it.  Also,  this  out  of  plane  mode  is  uncoupled  from  the 
planar  modes  (Ref  9)  and  hence  suffered  no  cross-feed  from  the  others. 

The  final  case  examined  was  that  of  calculating  the  gain  matrix  as 
a  function  of  time  based  upon  the  technique  of  Chapter  V.  The  singular 
value  decomposition  of  the  nonsquare  B  matrix  was  obtained  and  used  to 
get  (58).  Only  relocation  of  the  unstable  real  root  was  attempted. 

For  each  type  increment  in  the  period,  (63)  was  evaluated  to  check 
controllability  using  an  m  of  six  and  r  of  3.  The  values  from  (63) 
filled  a  three-by-six  matrix  which  was  inspected  for  nonzero  elements. 
For  each  increment  of  time,  one  of  the  three  columns  in  the  matrix  was 
nonzero.  The  controllability  condition  was  never  met  so  no  exact 
solution  was  obtainable.  This  was  verified  by  calculating  the  Poincare 
exponents  of  the  controlled  system  using  the  gain  from  (62).  Table  IV 
lists  the  values  obtained. 


TABLE  IV 


Poincare  Exponents  for  System  with 
Unmet  Control labil ity  Condition 

Desired  Location  Resultant  Location 


0.0 

+ 

0.077738i 

-0.000393 

+ 

0.077726i 

0.0 

- 

0.077738i 

-0.000393 

- 

0.077726i 

-0.1 

+ 

O.Oi 

0.19522 

+ 

O.Oi 

-0.17893 

+ 

O.Oi 

-0.17876 

+ 

O.Oi 

0.0 

+ 

0.084134i 

0.00011 

+ 

0.08231i 

0.0 

- 

0.084 134 i 

0.00011 

_ 

0.08231i 

from  this  table  it  is  seen  that  each  mode  was  effected.  The  fact 
that  the  unstable  real  root  maintained  its  positive  value  (it  actually 
became  more  positive)  verifies  the  uncontrollable  nature  of  the  system 
using  pole  placement.  No  case  was  developed  where  it  was  thought  the 
controllability  condition  would  meet. 

The  use  of  pole  placement  in  time  periodic  systems  is  highly 
effective  when  all  of  the  states  are  controllable.  Predictable  place¬ 
ment  of  poles  was  achieved  for  all  cases  in  a  completely  controllable 
system.  For  systems  with  uncontrollable  states,  such  as  positions  in 
the  restricted  four-body  problem,  predictable  placement  of  system  poles 


was  not  achieved. 
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